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Abstract 

Obesity is a major public health problem with a significant genetic component. Multiple DNA polymorphisms/genes have 
been shown to be strongly associated with obesity, typically in populations of European descent. The aim of this study was 
to verify the extent to which 6 confirmed obesity genes (FTO, CTNNBL1, ADRB2, LEPR, PPARG and UCP2 genes) could be 
replicated in 8 different samples (n = 1 1,161) and to explore whether the same genes contribute to obesity-susceptibility in 
populations of different ancestries (five Caucasian, one Chinese, one African-American and one Hispanic population). GWAS- 
based data sets with 1000 G imputed variants were tested for association with obesity phenotypes individually in each 
population, and subsequently combined in a meta-analysis. Multiple variants at the FTO locus showed significant 
associations with BMI, fat mass (FM) and percentage of body fat (PBF) in meta-analysis. The strongest association was 
detected at rs7185735 (P-value = 1.01 xlO" 7 for BMI, 1.80x10" 6 for FM, and 5.29x10" 4 for PBF). Variants at the CTNNBL1, 
LEPR and PPARG loci demonstrated nominal association with obesity phenotypes (meta-analysis P-values ranging from 
1.15x1 0~ 3 to 4.94 x 1 0~ 2 ). There was no evidence of association with variants at ADRB2 and UCP2 genes. When stratified by 
sex and ethnicity, FTO variants showed sex-specific and ethnic-specific effects on obesity traits. Thus, it is likely that FTO has 
an important role in the sex- and ethnic-specific risk of obesity. Our data confirmed the role of FTO, CTNNBL1, LEPR and 
PPARG in obesity predisposition. These findings enhanced our knowledge of genetic associations between these genes and 
obesity-related phenotypes, and provided further justification for pursuing functional studies of these genes in the 
pathophysiology of obesity. Sex and ethnic differences in genetic susceptibility across populations of diverse ancestries may 
contribute to a more targeted prevention and customized treatment of obesity. 
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Introduction 

Obesity, a state in which excess lipids accumulate in various 
body fat depots due to a chronic imbalance between energy intake 
and energy expenditure, is associated with many diseases such as 
type 2 diabetes mellitus, hypertension, coronary heart disease, and 
some cancers [1]. Although the true causes of the accelerating 
obesity epidemic have not been fully clarified, the prevalence of 
obesity continues to increase around the world. If current trends 
continue, it is estimated that by the year 2030 nearly 50% of adults 
in the United States will be clinically obese [2], and the world 
population will include 1.12 billion obese individuals [3] . Although 
the impact of environmental factors is likely to be significant, many 
studies have shown that body weight and obesity are strongly 
influenced by genetic factors, with heritability estimates often in 
excess of 50% particularly when derived from comparisons of 
identical and fraternal twins [4,5]. 

Most genetic studies of obesity have focused on body mass index 
(BMI) to evaluate whether a person is obese or not. However, BMI 
cannot distinguish fat mass from fat free mass [6,7]. Alternative 
measurements such as percentage of body fat (PBF) and fat mass 
(FM) are more homogeneous and reflect body fat content more 
accurately than BMI [8] . Few genetic studies of obesity, however, 
have utilized these latter measures of body fat content. 

In the past few years, extensive efforts had focused on the 
detection of obesity genes. In this regard, relatively few genes 
identified through preliminary linkage scans or candidate gene 
approaches have been confirmed to be truly associated with 
obesity by replication studies and other methodologies [9] . Among 
potential reasons for the failure to replicate most candidate genes 
or linkage peaks, one can cite small sample sizes and limited 
numbers of DNA variants upon which these studies were based. 
More recendy, under the common variant-common disease 
hypothesis, several genome-wide association studies (GWAS) and 
large-scale meta-analyses of multiple GWAS on obesity (mostiy 
using BMI as phenotype) have been conducted. Thus far, a total of 
58 genetic loci, all with small effect sizes, were found to be robusdy 
associated with obesity-related traits in multiple populations ([10]; 
see also: http://www.ncbi.nlm.nih.gov/gap/phegeni/). Among 
these genes, fat mass- and obesity-associated gene (FTO) stands 
out as the gene with the strongest significant association with 
obesity and it has been found to be associated with obesity in 
virtually all populations in which replication was attempted. 
Nevertheless, all genomic markers identified along with their 
putative genes have only been shown to have very small effects on 
BMI or the risk of obesity. Cumulatively, these genetic loci 
identified through GWAS account for less than 5% of the total 
heritability of BMI [1 1], leaving the vast majority of heritability yet 
unidentified. Presumably, additional variants/loci will eventually 
be detected with larger sample sizes combined with incorporation 
of rare variants, copy number variation markers, and other 
genomic and epigenomic features. 

The leptin receptor (LEPR), a single-transmembrane-domain 
receptor of the cytokine receptor family [12], acts with an 
adipocyte-specific hormone leptin that regulates adipose-tissue 
mass through hypothalamic effects on satiety and energy 
metabolism. Peroxisome proliferator-activated receptor gamma(P- 
PARG) is a transcription factor expressed abundantly in adipose 
tissue which involved in adipogenesis by activating adipocyte 
differentiation and mediating the expression of fat cell-specific 
genes [13]. Adrenoceptor beta 2, surface (ADRB2) is a major 
lipolytic receptor in human fat cells which plays a key role in 
regulating energy balance through both thermogenesis and lipid 
mobilization from adipose tissues [14]. Uncoupling protein 2 



(UCP2) is an inner mitochondrial membrane transporter which 
dissipates the proton gradient of inner mitochondrial membranes 
and releases stored energy as heat, thus has an important role in 
energy expenditure [15]. FTO gene encoding a nucleic acid 
demethylase plays a role in controlling feeding behavior and 
energy expenditure [16]. Catenin (cadherin-associated protein), P- 
like 1 (CTNNBL1) encodes a protein homolog to P-catenin 
responsible for cell-to-cell adhesion and Wnt-signalling [17]. 
Replication of previous findings in well-designed and statistically 
powered studies and fine-mapping the causal variants are essential 
to elucidate the importance of these six genes on obesity. The 
present study represents a meta-analysis of samples with diverse 
ancestries in which we attempt to replicate associations between 
sequence variants in six important candidate genes of obesity and 
body fatness phenotypes and to explore whether the same genes 
contribute to obesity-susceptibility in populations of different 
ancestries in large study samples. 

Materials and Methods 

Study Populations 

We utilized eight GWAS, four of which were "in-house" 
studies: (1) Quebec Family Study (QFS, n = 875, Caucasian 
Ancestry), (2) Omaha Osteoporosis Study (OOS; n = 998, 
Caucasian Ancestry), (3) Kansas-City Osteoporosis Study (KCOS; 
n = 2,283, Caucasian Ancestry) and (4) China Osteoporosis Study 
(COS; n= 1,624, Han Chinese Ancestry). Three studies(OOS, 
KCOS and COS) were originally designed to identify potential 
genes underlying osteoporosis. The other four were "external" 
studies deposited into the Database on Genotypes and Phenotypes 
(dbGaP) at the National Library of Medicine (http:/ /www.ncbi. 
nlm.nih.gov/gap/): (1) Framingham Heart Study (FHS; n = 2,786, 
Caucasian Ancestry), a longitudinal and prospective cohort 
comprising over 16,000 Caucasian subjects spanning three 
generations. Based on the first two generations of the FHS 
families, we identified 2,786 subjects with both BMI and FM 
information for use in this study. (2) Indiana Fragility Study (IFS; 
n= 1,478, Caucasian Ancestry), a quantitative and cross-sectional 
cohort comprising premenopausal Caucasian sister pairs. (3) 
Women's Health Initiative (WHI) Observational Study [18] 
African Sub-study (WHI-AA; n = 709, African Ancestry), (4) 
WHI Observational Study Hispanic Sub-study (WHI-HIS; 
n = 408, Hispanic Ancestry). Details regarding these studies have 
been published previously [19-21]. All studies were approved by 
their respective institutional ethics review boards. 

Phenotype Measurements 

Several obesity-related phenotypes were measured. These 
include BMI, body composition (FM and PBF) measured by 
dual-energy X-ray absorptiometry (DXA) scanners (either Lunar 
Corp., Madison, WI, USA, or Hologic Inc., Bedford, MA, USA) 
following the manufacturer's protocols or by underwater weighing 
with corrections for pulmonary residual volume. Covariates, 
including sex, age, age 2 , weight, height, and scanner ID (in 
WHI-AA and WHI-HIS), were screened with a stepwise linear 
regression model. Raw measurements were adjusted for significant 
covariates. To correct for potential population stratification, 
principal components (PCs) were computed and the first five 
PCs (i.e., PC1-PC5, explained >95% variation in each popula- 
tion.) derived from genome-wide genotype data were also included 
as covariates. Residual scores for each phenotype were normalized 
by inverse quantile of the standard normal distribution to impose a 
normal distribution on phenotypes which were then subjected to 
further analysis. 
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Genotyping and Quality Control 

All eight cohorts were genotyped using high-throughput SNP 
genotyping arrays (Affymetrix Inc., Santa Clara, CA; or Illumina 
Inc., San Diego, CA, USA) following their respective manufac- 
turer's protocols. Implemented in PLINK (http:/ /pngu.mgh. 
harvard.edu/~purcell/plink/), quality control criteria included 
the following: missing data <5%, SNP call rate >95%, and 
Hardy-Weinberg equilibrium (HWE) P-value> 1 .0 x 10~ 5 . For two 
family-based studies (i.e., FHS and IFS), all genotypes with 
Mendelian inheritance errors were set to missing. Details 
regarding the genotyping platforms, quality control, and data 
cleaning measures have been described previously [19-21]. 

Genotype Imputation 

To combine data across different genotyping platforms and to 
achieve a higher genome coverage, extensive genotype imputation 
was performed. Briefly, haplotypes of individual GWAS were first 
phased by a Markov Chain Haplotyping algorithm (MACH) 
(http://genome.sph.umich.edu/wiki/MaCH) [22], and untyped 
genotypes were then imputed by Minimac (http://genome.sph. 
umich.edu/wiki/Minimac), based on phased haplotypes, using the 
freely available haplotype data of the 1000 Genomes Project (as of 
August, 2010) as reference panels. Reference samples included 283 
individuals of European ancestry, 193 individuals of Asian 
ancestry, and 174 individuals of African ancestry. Imputation 
was performed by comparing the respective panel with the closest 
ancestry. For each GWAS, genotypes for untyped SNPs were 
imputed based on relevant population's haplotype reference panel. 
SNPs with imputation quality score (as assessed by r2.hat by Minimac) 
greater than 0.3 were retained in at least two studies, and with 
minor allele frequency (MAF) >0.05 in at least one study, were 
included for subsequent analyses. Prior to genotype imputation, 
strand orientations were checked and inconsistencies were 
resolved. Imputation results are summarized as an 'allele dosage' 
defined as the expected number of copies of the coded allele at that 
SNP (i.e., a fractional value between 0 and 2) for each genotype. In 
total, 2,954 genotyped or imputed autosomal SNPs spanning the 
FTO(1,275 SNPs), ADRB2(U9 SNPs), CTMMBL1{483 SNPs), 
LEPR{550 SNPs), PPARG{435 SNPs) and UCP2 (92 SNPs) genes 
were analyzed. 

Association Tests 

Association was tested in each study between directly-typed or 
imputed SNPs and obesity phenotypes under an additive genetic 
model. For each study of unrelated subjects (i.e., OOS, KCOS, 
COS, WHI-AA, and WHI-HIS), association was examined by 
fitting a linear regression model using MACH2QTL (http:/ /www. 
sph.umich.edu/ csg/abecasis/MACH/download/), in which allele 
dosage was used as a predictor of phenotype. For family-based 
samples (i.e., QFS, FHS and IFS), a mixed linear model was used 
in which the effect of genetic relatedness within each pedigree was 
taken into consideration [23]. 

Meta-analysis 

Conventional meta-analysis of individual studies was performed 
with weights proportional to the square root of the sample size 
using METAL software (http:/ /www.sph.umich.edu/csg/ 
abecasis/metal/) [24], and Cochran's Q statistic and f were 
calculated as measures of between-study heterogeneity. Random- 
effect meta-analyses were performed particularly for SNPs with Qi 
statistic P-value<0.05 or / value>50%. Meta-analysis of effect 
size (regression slope) of the identified candidate SNP and the 



forest plot were performed using Review Manager (RevMan) 
(http://ims.cochrane.org/ revman). 

Regional association plots of the most significant SNPs were 
generated using LocusZoom [25]. Wright's F-statistics (Fgx) was 
calculated with the R package to assess genetic differences among 
populations using the differences in allele frequencies. We 
interpreted the resultant F S x values based on Wright's suggested 
qualitative guidelines of F ST values as follow: F ST = 0-0.05 
indicating little population differentiation, 0.05-0.15 indicating 
moderate differentiation, 0.15-0.25 indicating great differentia- 
tion, and >0.25 indicating very great differentiation. 

Trans-ethnic Meta-analysis 

Given the multiple ancestry groups included in the study, we 
further performed a trans-ethnic meta-analysis using a recent 
MANTRA (Meta-Analysis of Trans-ethnic Association studies) 
software [26] as an alternative and robust analytic approach that 
accounts for the trans-ethnic nature of the various cohorts. In 
contrast to the traditional meta-analysis, the trans-ethnic approach 
in MANTRA accommodates between-population heterogeneity of 
associated variants and their effect sizes by allowing for allelic 
effects to be most similar between the most closely related 
populations. MANTRA adopts a Bayesian framework and 
assumes that studies from closely related populations are more 
likely to share a common true effect size, and the true effect size is 
allowed to vary across different population clades. Default settings 
were used in MANTRA. Evidence in favor of association of the 
trait with the variant was assessed by means of a Bayes' factor (BF). 

Results 

Baseline characteristics of the subjects in each cohort are 
presented in Table 1. A total of 2,954 SNPs (77.2%-91.1% 
imputed SNPs) covering the six targeted genes were meta- 
analyzed. A total of 45 SNPs located in the FTO gene 
demonstrated significant associations with obesity phenotypes 
using a Bonfferoni corrected significance level of P< 1.68x10 5 . 
Linkage disequilibrium (LD) analysis revealed that these 45 SNPs 
were in almost complete LD (r 2 >0.80) and were located within the 
same LD block of approximately 50 kb (Figure 1). Another 125 
SNPs out of the 1,275 total SNPs spanning the FTO gene in the 
present study showed marginal associations with obesity pheno- 
types (meta-analysis P-values ranging from 1.71x10 5 to 
4.99 xlO -2 ) (Results not shown). 

Table 2 shows the top 15 FTO SNPs associated with obesity 
traits. For most variants, little heterogeneity among cohorts was 
observed, except for rs9922708 for which moderate heterogeneity 
was observed (I 2 = 43% and Qp- Vil | ue = 0.091). The direction of 
effects was consistent across studies except for three SNPs 
(rs9922708, rs 178 17449 and rs7206790). The trans-ethnic meta- 
analysis using MANTRA closely mirrored the results from the 
standard meta-analysis using METAL. The most significant 
association was observed for FTO rs7185735 (P-val- 
ue=1.01xl0" 7 for BMI, 1.80xl0" 6 for FM, and 5.29xl0" 4 
for PBF) (Table 2), and allele-specific OR (95% CI) reached 1.20 
(1.12-1.28) for BMI, 1.18(1.10-1.28) for FM and 1.13 (1.05-1.21) 
for PBF (Figure 2). The frequency of minor allele G at this 
imputed SNP ranged from 0.11 (Chinese) to 0.47 (Caucasian) in 
different ethnic groups. The forest plot for this SNP indicated that 
its strongest association was observed in QFS (OR= 1.51 for BMI, 
1.49 for FM and 1.37 for PBF). Carriers of the G allele had higher 
BMI, FM and PBF values (Figure 2). 

For the other five genes, 26 SNPs, 1 1 SNPs and 33 SNPs from 
CTMNBL1, PPARG and LEPR showed marginal significance of 
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it associations with P-values ranging between 

-g- 4.94xlO _2 (Results not shown). The strongest associations for the 

K CTNNBL1, PPARG and LEPR gene markers were: CTNNBL1 

| rs45500793 (P-value = 2.68x 10~ 3 for BMI, 6.34xl0" 3 for FM, 

| and 2.06 x 10" 2 for PBF), PPPP rs9436744 (P-value = 3. 16 x 10~ 3 

S. for BMI, 2.33 x 10" 2 for FM, and 7.55 x 10" 2 for PBF) and PPARG 

£ rsl0222537 (P-value = 1.50xl0~ 3 for BMI, 1.19xl0" 2 for FM, 

s and 6.77 xlO -2 for PBF). We did not detect association between 

i/,- any variants in the ADRB2 and UCP2 genes and obesity 
phenotypes. 

"g To identify potential sex-specific loci, we performed a series of 

" meta-analyses stratified by sex. In female specific samples, 25 FTO 

o SNPs were significantly associated with BMI and FM, all of which 

§- overlapped with the findings obtained with the combined sample. 

$ There was no evidence of significant associations in the male 

ra specific samples (Table 3). However, a nominally significant 

g association was revealed at rsl6952725 of FTO in males (P- 

g value = 6.87 xlO" 4 for BMI and P-value = 8.55 x 10" 3 for FM), 

y_ but not in the female specific sample (P-value = 0.55 for BMI and 

f" P-value = 0.43 for FM) (Table 3). We estimated the statistical 

£ power at various locus heritabilities for the total sample, males and 

o females (Figure 3). Our total sample and female sample have > 

g- 90% power to detect a variant explaining >0.40% of heritability. 

2 However, the male sample only have 30% power to detect a 
>, variant explaining >0.40% of heritability. 

V Finally, we examined potential differences of SNP association 

c with obesity phenotypes among the four ethnic groups (Chinese, 

. Caucasians, African Americans and Hispanic-Americans). In 

g Caucasians, we found that 35 FTO SNPs were significandy 

* % associated with BMI and FM (P< 1 .68 x 10~ 5 ); 28 of these 35 SNPs 

"§ 5. overlapped with those derived from the analysis performed on the 

j2 ^ combined samples. The most significant association was observed 

£ J for rsl7817288 in PTO (P-value= 1.96x 10" 7 for BMI, 2.76x 10" 5 

| :| for FM, and 2.39 x 10~ 3 for PBF). However, these significant FTO 

K j= SNPs in Caucasians showed no or weak evidence of association in 

0 ±± 

™ S non-Caucasian populations (Table 4, Figure 1). Therefore, the 

1 y association results of FTO obtained with the combined samples 
o c . . 

| were mostly contributed by the Caucasian group. For example, P- 

8 | value was 1.01 xlO" 7 for the most significant SNP rs7 185735 in 

§sd the combined samples for BMI, 4.90x10 7 in Caucasians, 

£ j 7.47x10 2 in Chinese, 6.05x10 1 in African-Americans and 

= 5 4.45x10 1 in Hispanic-Americans. In non-Caucasian popula- 

£ § tions, we only found marginally significant associations. Two FTO 
SNPs (rs7201444 and rs 1 3335 146) specific in Hispanic -Americans, 
14 FTO SNPs and one LEPR SNP(rs9436299) specific in African 
Americans and one CTNNBL1 SNP (rs45500793) specific in 

a < Chinese showed marginally significant associations (P-value< 

|.| lxlO" 3 ) (Table 4, Figure 1). Table 4 listed 2 SNPs specific in 

S ■£ Hispanic-Americans, 15 SNPs specific in African Americans, 1 

S! ~ SNP specific in Chinese and top 15 SNPs specific in Caucasians. 

^ "| To assess genetic differences for these ethnic-specific SNPs (listed 

3 ^ in Table 4) among the eight populations, F S t were calculated. 
Nine FTO variants (rsl861869, rs7186521, rsl7817288, 
rs8044769, rsl 1075987, rs9935401, rs8051591, rs7193144, 



a c 



— '-I 

I «e | rs8043757) and CTNNBL1 SNP (rs45500793) showed little 

population differentiation (F ST = 0-0.05). F ST of other SNPs 
J 5 J ranged from 0.05 to 0.31, suggesting moderate/great genetic 

differentiation. FTO rs7201444 showed very great differentiation 
j! tt E with the frequency of minor allele A ranging from 0.001 (Chinese) 

3 1= ^. to 0.365 (African-Americans). 

S 2 ^ In addition, we examined the association of previously reported 

i ? 5 SNPs of the six genes among our dataset (Table 5). All of the 10 

Hi previously reported FTO SNPs (rsl421085, rsl558902, 

rsl7817449, rs9941349, rs8050136, rsl558902, rsl 121980, 
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Figure 1. Regional plots of /TO gene in African (A), Chinese (B), Caucasian (C) and Hispanic populations (D), respectively. SNPs are 
plotted by position on the chromosome against association with BMI (— Iog10 P-value). Recombination rates (from HapMap) are plotted in blue to 
reflect the local LD structure. The SNPs surrounding the most significant SNP (in purple) are color coded to reflect their LD with this SNP (r 2 values 
from the 1000 Genomes Mar 2012 AFR, ASN, EUR and AMR data, respectively). 
doi:10.1371/journal.pone.0096149.g001 



rs7202116, rs9939609 and rs9930506) showed significant/nomi- 
nal associations with BMI, FM and/or PBF. Two PPARG SNPs 
(rs 180 1282 and rs3856806) and 3 CTNNBL1 variants (rs60 13029, 
rs 16986921 and rs6020712) showed nominal associations with 
obesity traits. 

Discussion 

In this study we investigated the associations between 6 genes 
and obesity phenotypes, including direct measures of body fatness, 
using a dense set of variants genotyped or imputed in 11,161 
subjects from four different ethnic backgrounds (Chinese, Cauca- 
sian, African-Americans and Hispanic-Americans). Our findings in 
this multiethnic population confirmed the importance of the FTO 
gene for obesity risk in humans. A total of 45 SNPs located in the 
FTO gene showed significant associations with the obesity 
phenotypes. The FTO protein affects demethylation of nuclear 
RNA in vitro [27], but whether the efficiency of this process 
depends on the FTO genotype or how this may be related to the 
observed effects on BMI or body fatness is not clear. 

At present, the strongest associations between FTO SNPs and 
BMI belong to intronic SNPs, which might have a role in the 
regulation of FTO and/ or nearby genes. It is critical to recognize, 
however, that associated SNPs are not necessarily causal SNPs 



underlying the association, and that the functional variants are still 
unknown. A recent study reported that rs72021 16 G allele creates 
a CpG site along with other variants in perfect linkage 
disequilibrium with it [28], and these risk alleles may potentially 
have increased DNA methylation. Bioinformatic analyses also 
revealed that six FTO SNPs (rs 1 1 6420 15, rs 1 78 1 7497, rs375 1 8 1 2, 
rsl7817964, rs62033408, and rsl421085) were located within 
candidate intronic regulatory elements and that two SNPs 
(rsl 1642015 and rsl421085) were predicted to have allele-specific 
binding affinities for different transcription factors [29]. Specifi- 
cally, the T allele at rsl 1642015 binds Paired box protein 5 
(PAX5) while the C allele at rsl 42 1085 is predicted to have a 
substantially reduced binding affinity for Cut-like homeobox 1 
(CUX1). Further investigation is warranted to identify potentially 
functional SNPs and the mechanisms by which various alleles at 
FTO influence the level of adiposity. 

We detected moderate associations for multiple variants from 
the CTNNBL1, PPARG and LEPR genes, but did not detect 
association with any variants in ADRB2 and UCP2. Our data 
confirmed the associations with adiposity for three previously 
reported CTNNBL1 variants (rs60 13029, rsl 6986921 and 
rs6020712) [17]. This finding is consistent with the results from 
Andreasen and co-workers, who in a study comprising 18,014 
Danish participants found that the CTJVJVBL1 rs60 13029 T-allele 
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Figure 2. Forest plot of the association of rs71 85735 and BMI (upper), FM (middle) and PBF (lower). 
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Figure 3. Power estimation for the total sample, females and males. 

doi:1 0.1 371 /journal.pone.00961 49.g003 



and the rs6020846 G-allele confer an increased risk of developing 
obesity, especially morbid obesity [30]. Interestingly, a German 
population-based study (KORA) [31] failed to demonstrate 
association of CTMMBL1 variant rs60 13029 with obesity or BMI. 
The discrepancy might be due to population-specific differences 
and/or insufficient statistical power. For the most frequendy 
studied coding variants of UCP2 [32], ADRB2 [33], PPARG [34] 
and LEPR [35] genes, only rs 180 1282 and rs3856806 in PPARG 
gene showed nominal associations with BMI and/ or PBF. 

There are substantial ethnic differences in the prevalence of 
excess body weight and obesity. Although differences in lifestyle 
are likely to account for some of the observed differences, genetic 
variability could also play a major role. Among the four ethnic 
groups (Chinese, Caucasian, African-Americans and Hispanic 
populations), we found that 35 SNPs in intron 1 of FTO gene were 
significandy associated with indicators of obesity in Caucasians. All 
of the SNPs overlapped with those reported in a previous meta- 
analysis study [36]. There was very limited or no evidence for 
associations between these SNPs and adiposity in the other ethnic 
groups. On the other hand, two FTO SNPs (rs7201444 and 
rsl 3335 146) in intron 8 specific in Hispanic Americans, 11 FTO 
SNPs (rs8043738, rs8063472, rs8058460, rsl077129, rs8048396, 
rsl 86 1869, rsl 72 17 144, rs2892469, rsl 86 1868, rs7 186521, 
rs7 184874) in intron 1 specific in African Americans showed 
moderate evidence for associations that were not reported by 
previous studies. While several previous studies reported associa- 
tion between FTO SNPs and obesity-related phenotypes in 
Hispanic Americans [37-39], African Americans [29,40] and 
Asian populations [41 — 44], these studies were relatively small in 
sample sizes and showed mixed results. Wing et al. [39] genotyped 



26 SNPs in intron 1 of FTO in 373 Hispanic Americans and 
observed associations between BMI and several SNPs that were 
previously reported to be associated with obesity (rs9939609, 
rs8050136, rsl 121980, rsl421085, rsl7817449 and rs3751812), 
and four other SNPs (rs8047395, rsl0852521, rs8057044 and 
rs8044769). Song et al. [37] and Scuteri et al. [38] replicated 
associations of BMI with the well-known rs9939609 [37] and 
rs9930506 [38] in Hispanic Americans. However, their findings 
were not confirmed in our Hispanic Americans. The two 
Hispanic-American specific SNPs found in our study were in 
intron 8 of FTO. Interestingly, a detailed mapping (262 tag SNPs 
across the entire FTO gene) in individuals of African American 
descent demonstrated significant association for seven SNPs 
(rs708262, rsl 1076017, rsl6952725, rs9932411, rs7191513, 
rs2689269, rsl 6952987) in intron 8 [40]. Our study together with 
Adeyemo et al. suggests the intron 8 of FTO is a second site, in 
addition to intron 1 , playing a role in the association between FTO 
and indicators of obesity. 

Previous studies in African Americans showed either very 
limited or no evidence for associations with the SNPs initially 
reported in European populations, such as rs9939609, rsl 121980, 
rsl 78 17449 and rs8050136 [29]. More comprehensive evaluations 
of the FTO locus identified several other variants, rs56137030 
[29], rs3751812 [45,46], rsl 108102 [39] and rs8057044 [47] that 
showed associations with BMI in African Americans. Among these 
variants, only rsl 108102 (P- value = 6.15 x 10 ) replicated in our 
African American sample. Among the 14 FTO SNPs specific to 
African Americans, rsl 86 1868 was previously found to be the 
strongest SNP associated with BMI in Old Order Amish 
individuals with low physical activity (P<0.001) [48]. rsl 86 1869 
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and rs7 186521 were found to be associated with weight and waist 
circumference in 843 unrelated individuals from an island 
population in the eastern Adriatic coast of Croatia [49]. The 
allele frequencies for most of the SNPs indicated moderate 
population differentiation, and thus could potentially lead to 
variable genetic impact on obesity phenotypes across populations. 
In summary, our study replicated the associations of FTO intron 1 
variants with BMI, FM and PBF in Caucasians, and confirmed the 
importance of variants in intron 8 as well. Therefore, our results 
suggest that there are ethnic differences with regard to the effects 
of FTO on obesity and body fatness. 

Marked sex differences in the prevalence of obesity between 
women and men suggest sex-specific genetic impacts on obesity 
risk, even though lifestyle and dietary factors are also likely to 
contribute to these differences. Furthermore, genetic factors may 
interact with levels of physical activity to modify obesity risk [50] . 
To determine if the effect of genetic variants differed by sex, we 
stratified by sex and found that 25 FTO SNPs were specifically 
associated with obesity in females only. This is consistent with a 
few studies that have shown sex-differences in the heritability of 
BMI and fat percentage [51,52]. In contrast, a nominally 
significant association was revealed at rsl6952725 of the FTO 
gene in males only, but not in females. Our findings are further 
supported by a recent study in children which found that the FTO 
variant rs9939609 showed association with obesity and BMI 
among girls but not among boys [53]. A recent genome-wide 
association study modeling the effect of genotype-by-sex interac- 
tion on obesity phenotypes demonstrated sex-influenced associa- 
tions between genetic variation at the LYPLAL1 locus and obesity- 
related traits [54]. This study suggests that FTO may be a gene 
playing a role in the commonly observed sex-dimorphism in 
adiposity. 

The majority of obesity loci have been discovered through 
GWAS in individuals of European descent, and, more recendy, in 
Asians as well. Okada et al. performed a GWAS with BMI in 
62,245 East Asian subjects and observed a significant association 
with FTO rsl2149832 (P-value = 4.8 x 10~ 22 ) [55]. However, our 
study failed to replicate the association of FTO rsl2149832. 
Several other studies replicated FTO rs9939609 [41-43,56,57] and 
rs8050136 [44] with obesity and BMI in Chinese or Asians. Our 
study showed weak evidence of association with FM and PBF for 
rs9939609 and rs8050136 (P-values range from 9.53 xl0~ 3 to 
3.23 xlO -2 ), but not with BMI. Lu and Loos reviewed the 
transferability of 36 GWAS identified BMI-associated SNPs 
between European and East Asian ancestry populations using 
both SNP-to-SNP and locus-wide comparisons [32]. SNP-to-SNP 
comparisons suggest that of the 32 SNPs found in European 
populations, 8 SNPs were non-polymorphic and another 6 showed 
no convincing evidence of association in East Asians and two of 
the four loci identified in East Asians showed some evidence of 
transferability to European populations [32]. However, locus-wide 
analyses suggested the more extensive transferability. For example, 
CDKAL1 rs9356744 showed genome-wide significant association 
with BMI in East Asians but not in European populations (P- 
value = 0.19). When examining the locus surrounding rs9356744, 
however, other SNPs that are not in LD with rs9356744 were 
actually associated with BMI in European populations (P-values< 
10 -3 ). Replication studies in other populations that scan the entire 
gene (rather than just one or a few SNPs) may lead to the discovery 
of other important genetic variants. Studying populations of 
different ancestries (especially those with smaller LD) could help 
fine-map disease or trait loci, eventually pinpointing the causal 
variant (s) and gene(s). Our study used a comprehensive approach 
by thoroughly examining six genes with dense SNP coverage to 



search for population-specific and/or shared obesity loci in 
multiple ancestry groups. 

An important limitation of the current study is the relatively 
small numbers of male subjects (29% of total subjects) and subjects 
from the three non-Caucasian ethnic groups (Chinese, African and 
Hispanic populations, 25% of the total subjects). These low 
numbers may have contributed to the nominal associations 
detected in these sub-groups. Nevertheless, the present study 
provided additional evidence supporting the presence of ethnic- 
and sex-differences for some prominent obesity variants and 
perhaps genes. 

In conclusion, we have attempted to replicate previously 
reported associations between multiple DNA common variants 
relating to 6 obesity genes in populations representing four 
ethnicities. We were able to find confirmatory evidence for 
contributions of FTO, CTMNBL1, LEPR and PPARG related 
genomic variants to human variation in adiposity. In particular, 
FTO variants showed sex-specific and ethnic-specific associations 
with adiposity traits. 
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